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ABSTRACT 



Context. Gamma-ray bursts (GRBs) and their early afterglows ionise their circumburst material. Only high-energy spectroscopy 
therefore, allows examination of the matter close to the burst itself. Soft X-ray absorption allows an estimate to be made of the total 
column density in metals. The detection of the X-ray afterglow can also be used to place a limit on the total gas column along the line 
of sight based on the Compton scattering opacity. Such a limit would enable, for the first time, the determination of lower limits on 
the metallicity in the circumburst environments of GRBs. 

Aims. In this paper, we determine the limits that can be placed on the total gas column density in the vicinities of GRBs based on the 
Compton scattering. 

Methods. We simulate the effects of Compton scattering on a coUimated beam of high energy photons passing through a shell of 
high column density material to determine the expected lightcurves, luminosities, and spectra. We compare these predictions to 
observations, and determine what limits can realistically be placed on the total gas column density. 

Results. The smearing out of pulses in the lightcurve from Compton scattering is not likely to be observable, and its absence does not 
place strong constraints on the Compton depth for GRBs. However, the distribution of observed luminosities of bursts allows us to 
place statistical, model-dependent limits that are typically < 10^' cm"^ for less luminous bursts, and as low as - 
luminous. Using the shape of the high-energy broadband spectrum, however, in some favourable cases, limits as low as - 
can placed on individual bursts, implying metallicity lower limits from X- and gamma-rays alone from up to 0.01 Z/Zq. At extremely 
high redshifts, this limit would be at least 0.02 Z/Z©, enough to discriminate population 111 from non-primordial GRBs. 



Key words. Gamma-ray burst: general - Gamma-ray burst: individual: GRB 050904 - early Universe - 
stars - Galaxies: ISM - Stars: Population III 
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1. Introduction 

A major issue in modern cosmology is our understanding of 
the metal-enrichment history of the universe: how, when and 
where the non-primordial elements were synthesised. Most av- 
enues open to us to investigate this question at high redshifts rely 
on bright sources such as active galactic nuclei and gamma-ray 
bursts (GRBs) as back-lighting for absorption spectroscopy of 
distant galaxies. Gamma-ray burst afterglows are bright enough 
to examine the contents of star-forming galaxies by absorption 
spectroscopy. Most detailed information so far has come from 
optical and ultraviolet (UV) observations. Due to the strong ion- 
ising effect of the burst and early afterglow on the matter sur- 
rounding a GRB however only high-energy spectroscopy can 
probe the matter close to the burst. X-ray spectra allow exam- 
ination of th is matter, which we now know lies within only a 
few par secs jWatson et al.ll2007^_ and is ionised to a fairly high 
degree d^ochaska et al. 2008; Schadv et al.ll2010l) and is there- 
fore invisible at optical and UV wavelengths. 

Photoelectric absorption by inner shell electrons absorb X- 
rays and this allows an estimate to be made of the total column 
density in metals. The total gas column is unavailable, however, 
since the hydrogen, which is ionised, is transparent to the X-rays 
even at moderate column densities. However, at extreme column 
densities, ionised hydrogen is no longer transparent to high en- 



ergy photons since X-ray scattering off the free electrons has 
a cross section approximated by the Thomson cross section at 
low energies, <tt ~ 6.65 x 10"^^ cm^. This fact has been used 
to derive limits on the column density of ionis ed hydrogen in 
the im mediate surroundings of the GRB 050904 (ICampana et al.l 
|2007|) . and to exclude Comptonisation as the origin of very h ard 
gamma-ray emission in some GRBs ('Ghirlan da et al.l 20031) . It 
has also recently been suggested by[Campana e t al.l ( 2010al) that 
the detection of the X-ray and gamma-ray emission from GRBs 
could be used to place a limit on the total gas column along the 
line of sight based on the fact that at very high column densi- 
ties Compton scattering will essentially eliminate all of the X- 
and gamma-ray emission and that this could be useful for a fu- 
ture GRB mission. Such a limit would allow us to estimate lower 
limits to the metallicity in the immediate environments of GRBs 
by providing an upper limit to the total gas column density. At 
high redshifts, not only could this allow us to determine metal- 
licities without the problems associated with UV spectra, but it 
would also allow us to determine the metal abundances right in 
the hearts of star-forming regions. 

Future X-ray missions equipped with large effective area, 
high resolution spectrographs, would allow redshifts to be de- 
termined directly from X-rays without reference to optical or 
UV light, enabling far more complete and unbiased estimates of 
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GRB environments to be obtained. However, without this tech- 
nique, they could never determine the metalHcities of the re- 
gions they probe. In pl acing this Compton-thick limit, however 
ICampana et al.l (l2010ah assume for simplicity that this limiting 
column density is approximately one optical depth t = I/ctt. 
But given that the observed column densitie s of GRBs span a 
range fr om < 10^' cm ~^ up to ~ lO^^'cm"^ (lEvans et al 
FCampana et al.ll2010bb . the value of the limit that can be placed 
on the total column density makes a substantial difference to 
how interesting the metallicity lower limit obtained will be. In 
this paper we attempt to quantify, via Monte Carlo simulations, 
the limits that can placed on the total column density and hence 
metallicities of the circumburst medium from the Compton- 
transparency of GRBs and their afterglows. 

2. Simulations 

2.1. Burst modeling 

The radiative transfer (RT) calculations are conducted using a 
modified version of the Monte Carlo RT code MoCaLaTA, orig- 
inally designed for simulating the scattering of Lye? photons 
in the interstellar and intergalactic media (Laursen et al. 2009). 
Although the code is capable of assuming an arbitrary distribu- 
tion of gas density, temperature, velocity field, etc., taken e.g. 
from a hydrodynamic simulation, for the purpose of the present 
calculations the burst is modeled as a central source emitting 
a power-law spectrum of photons with energies E between 0. 1 
keV and 1 MeV and a photon spectral index of F = 2. That is, 
the emitted spectrum is given by 

Fe„(£) = A£-^ (1) 

where A is a constant. The source is surrounded by a thick, spher- 
ical shell of electron column density Ne, with a temperature T 
and velocity field Vbuik- Photons are emitted in two narrow cones 
of opening angle § = 5°, taken to lie along the z axis. 
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Fig. 1. Lightcurves for a range of electron column densities A^^,. 
Only the prompt emission will be detectable with current instru- 
mentation for cosmological GRBs. The inset is a zoom-in of the 
prompt emission. The vertical dashed line marks the maximum 
time lag for singly scattered photons, as given by Equation|4] 



MoCaLaTA has already been tested extensively against vari- 
ous analytical solutions; the additional subroutines implemented 
in this study were tested correspondingly. Furthermore, a variety 
of temperatures, velocity fields, density fields, and opening an- 
gles were tested. However, as will be shown in Section Sect. l3.1l 
observationally only the prompt emission is of interest, which 
depends almost exclusively on the column of intervening gas and 
not on its actual structure 



2.2. Radiative transfer 



3. Results 



While at low energies a photon interacting with an electron is 
scattered in a more or less random direction, at high energies 
the phase function is characterised by a significant probability 
of forward scattering. Moreover, the cross section crc decreases 
with increasing photon e nergy. These effects ar e given by the 
differential cross section (iKlein & Nishinairi929h 



— ^ = -r;P^(E, 0) P(E, 0) + 1 + cos" 

dQ. 2 ' ^ '\ ^ P(E,0) 

where r^, is the classical electron radius, and 



P(E, 0) 



1 



1 - q;(1 - cos 0) 



(2) 



(3) 



is the energy ratio of the incident and the scattered photon, with 
a = E/nieC^ being the photon energy in terms of the electron rest 
energy. 

The photons are followed as they scatter stochastically 
through the medium (or escape freely). Equations |2] and |3] are 
expressed in the reference frame of the electron. Thus, at each 
scattering the energy of a photon is first Lorentz transformed to 
the reference frame of the electron, given by the sum of Vbuik and 
a thermal velocity drawn from a Gaussian distribution. At high 
energies the change in energy is dominated by the Compton ef- 
fect described by Eq. [3j however, at lower energies the Doppler 
shift induced by high velocity electrons may become important. 



3.1. Liglitcurves 

Sampling the photons according to the time it takes to reach the 
observer, a lightcurve can be obtained, giving the photon flux as 
a function of time. Fig. [T] shows such lightcurves for a range of 
electron column densities. The curves are characterised by three 
parts of different physical origin: The initial peak — the "prompt" 
emission — consists of photons that escape the medium directly, 
without scattering, and hence its height is simply specified by the 
optical depth along the line of sight to the burst. The following 
intermediate phase consists of photons that are scattered once, 
and thus its extension is given by the size of the cloud and the 
opening angle of the jet; for a spherical cloud of radius R = 
1.25pc and opening angle i? = 5°, this phase lasts 



At = -(1 - cost?) ^ s^ie^ 

c 



(4) 



Finally, there is an afterglow of photons that scatter several times 
and thus perform a random walk out of the cloud. Since the num- 
ber of steps in a random walk scales with N^, but the step size 
decreases linearly with A^^ for a fixed physical size of the cloud, 
the maximum of the afterglow approaches a time proportional to 
N,. 

For high values of column density, the above considera- 
tions become somewhat inaccurate, since at each scattering, a 
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Fig. 2. Fraction of emitted photons that are transmitted directly 
through a cloud of column density Ne for the full energy range 
(red solid), low-energy photons [blue dashed), and high-energy 
photons {green dotted). 

photon — especially a high-energy photon — loses energy, thus 
resulting in a slowly decreasing electron cross section. 

Comparing the magnitude of the observed simulated prompt 
peak with the energy emitted by the burst, constraints can be 
put on the Compton thickness of the circumburst material. For 
the emitted spectrum in Equatior|TJ the observed spectrum of the 
prompt emission will be 



Fobs(E)^AE-' e 



(5) 



The fraction / of the flux observed as prompt emission is 



thus 



jj;"£-re-'v<-c(£)j£ 
J^-^E-^dE 



(6) 



This fraction is shown as a function of column density in Fig|2] 



3.2. Spectra 

Due to the decreasing Compton cross section with photon en- 
ergy, an emitted spectrum of constant F will have its high-energy 
end transmitted more efficiently than its low-energy end, re- 
sulting in a harder slope at high energies (Fig. [3]). In practice, 
this means that the spectral slope will change as a function 
of energy. Fig. |4] shows how the observed spectral slope Fobs 
changes with optical depth. Assuming an intrinsically energy- 
independent spectral slope F, if a broad energy range is observed 
the change in Fobs as a function of energy reveals the optical 
depth of the intervening cloud. For example, for F = 2 a col- 
umn density of A^^. ^ 10^^ cm"^ will result in an observed slope 
Fobs ^ 1 at high energies. 



4. Discussion 

To determine limits on the Compton depth of a given burst we 
have three diagnostics at our disposal: the lightcurve, the to- 
tal luminosity, and the spectrum. We have shown that the flux 
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Fig. 3. The model input and output spectra for the simulated 
GRB afterglow for a range of foreground column densities. 




Fig. 4. Change AF in the observed spectral slope Fobs as a func- 
tion of column density A^^., at various energies. 



in the lightcurve (and the fluence over the observed timescales 
of GRBs), is dominated by the directly transmitted component 
even up to high column densities (> 3 x 10^^ cm"^, see Fig.[TJ. 
This means that we cannot use the smearing of the shape of 
the lightcurve to discern even moderately high column density 
bursts. 

The limits that can be placed on the column density based 
on the luminosity of a burst are not very precise because we do 
not know intrinsically what the apparent luminosity of a burst 
was. We do, however, have a distribution of equivalent isotropic 
luminosities of GRBs, and we know that burst s with luminosi- 
ties a bove ~ 2 x 10^'* erg are rarely detected (iGhirlanda et al.l 
I2008h . If we assume that there is no bias relating luminous bursts 
with high column density sightlines, it suggests that any burst is 
likely to be intrinsically less luminous than this value. We can 
therefore use the luminosity of a burst to determine the diminu- 
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tion in flux possibly attributable to Compton scattering. In Fig.|2] 
we show the fraction of emitted photons in the prompt emis- 
sion as a function of column density. In the general case at high 
redshift, we are unlikely to detect any burst more than two or- 
ders of magnitude fainter than 10^"* erg with ^wi/f-BAT, and so 
the column density for high-redshift Swift bursts must always 
be < 10"^ cm"^ (see Fig.|2]i. In the specific case of the z - 6.3 
GRB 050904, this would imply a maximum column density of 
only a few times 10^'' cm"^ because of the high apparent lumi- 
nosity of this burst. However, individual cases are always vul- 
nerable to the argument that they may have an exceptional in- 
trinsic luminosity. Even in the general case, it is possible that 
z > 6 GRBs may be inherently more luminous than lower red- 
shift GRBs, weakening this argument further It might be possi- 
ble to predict the apparent luminosity of a burst based on a num- 
ber of correlations of burst properties with luminosity, the most 
well-regarded of which is the peak e nergy-isotropic lumino sity 
correlation (the "Amati relation", see iGhirlanda et al.ll2b08l and 
references therein). However, at the moment none of these cor- 
relations is very tight and their general validity is controversial. 
We therefore turn to the spectra to examine whether they can 
provide more robust constraints. 

The spectra show a clear deviation from a power-law at 
column densities above ~ 10^^ cm"^, if we have simultaneous 
broadband coverage from a few keV to a few hundred keV. 
Currently, we do not typically obtain simultaneous soft and hard 
X-ray detections of a GRB, though there are a few exceptions 
where the burst was long and the Swift response time short. A 
good example is, in fact, GRB 050904 where the spectrum of 
the late pro mpt phase is detected in t he range 2-700 k eV in the 
rest frame dWatson et al. 2006; Cus umano et al ll2007h . 

In general, at 10^^ cm"^, the spectral slope has changed 
by 1 in the 0.1-1 MeV range (Fig. |4]l. Such a strong change 
in the slope of a burst would be readily discernible in late 
prompt data with simultaneous X- and gamma-ray spectra, and 
it seems likely that this spectral characteristic represents the 
strongest constraint on the column densities. It should be cau- 
tioned that in most GRBs strong intrinsic spectral evolution of 
the burst occurs, causing the low energy data to change spec- 
tral slope substantially. In the case of GRB 050904, the slope 
changes in the first few minutes from F ~ 1.2 to F ~ 1.9. 
Such early spectral evolution is common in GRBs. It is there- 
fore essential that in looking for this effect the high and low 
energy data be simultaneous. Using GRB 050904, we note that 
the BAT and XRT contemporaneous spectra have compatible 
power-law s lopes, with uncertaint ies of ~ 0.1 on their spec- 
tral indices (ICusumano et al.ll2007h . Disregarding the relatively 
small flux and cross-calibration uncertainties, this similarity of 
the observed spectral slopes between the XRT and BAT yields a 
3cr limit on the column density of ~ 5 x 10~^ cm"^. This is a tight 
con straint, but still three t imes as high as the t - \ limit imposed 
by (ICampana et al.ll2007l) and yields therefore a metallicity limit 
of > 1% of the solar value (following their analysis of the X- 
ray absorption). It is worth noting that the d etection of the high 
column density in GRB 050904 is disputed dButler & Kocevskil 
l2007h . and there is the risk of confusing intrinsic spectral cur- 
vature for photoelectric absorption where we do not have high 
signal-to-noise ratio data at multiple epochs. 



primordial GRB would be useful. At z = 9-13, for exam- 
ple, the lowest feasible detection of photoelectric absorption is 
~ 10^^ cm"- in units of equivalent hydrogen column density at 
solar metallicity. Assuming a similar limit on the total electron 
column density as found for GRB 050904, of ~ 5 x 10^"* cm"^, 
where soft X-ray absorption is detected, this would imply a 
metallicity limit at least > 0.02 Z/Zq in the ~ 1 pc environment 
of the GRB. This should be sufficient to exclude a population III 
star as the GRB progenitor 

5. Conclusions 

We have examined the prospect of using only the high-energy 
emission of GRBs to place a lower limit on metallicities around 
them; the soft X-ray photoelectric absorption providing a mea- 
sure of the total metal column density and the Compton thick- 
ness limit a maximum column density of electrons and hence 
hydrogen. We find that smearing of the lightcurves does not pro- 
vide strong constraints on the column densities, while apparent 
luminosity and deviations from a power-law spectral shape pro- 
vide stronger constraints. The spectral constraints are more reli- 
able, but require the assumption of a single power-law spectral 
shape and simultaneous gamma- and X-ray coverage. The ap- 
parent luminosity constraints are more readily applicable, but 
are less certain, and rely heavily on an assumption that we know 
the distribution of apparent luminosities of GRBs, which is espe- 
cially uncertain at high redshift. We find for the individual case 
of GRB 050904, the best limit that can be obtained from spectra 
corresponds to ~ 5 x 10^'^cm"^, three times what was previ- 
ously assumed. While the results are not very constraining for 
most GRBs, with bright, long-lasting bursts, with detections of 
very high metal column densities, the method could be used to 
exclude a population III progenitor for a high redshift burst. 
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4.1. Population III progenitors 

At very high redshifts, GRBs may be formed from primordial 
stars - population III stars. If a GRB occurs at very high red- 
shift, a discriminator between a population III GRB and a non- 



